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InPBi was predicted to be the most robust infrared optoelectronic material but also the most difficult to 
synthesize within In-VBi (V = P, As and Sb) 25 years ago. We report the first successful growth of InPBi 
single crystals with Bi concentration far beyond the doping level by gas source molecular beam epitaxy. The 
InPBi thin films reveal excellent surface, structural and optical qualities making it a promising new III-V 
compound family member for heterostructures. The Bi concentration is found to be 2.4 ± 0.4% with 94 ± 
5% Bi atoms at substitutional sites. Optical absorption indicates a band gap of 1.23 eV at room temperature 
while photoluminescence shows unexpectedly strong and broad light emission at 1.4-2.7 \im which can't be 
explained by the existing theory. 

Compound semiconductors play a pivotal role in modern optoelectronics whose energy gap and lattice 
constant can be tailored to meet various device requirements. While semiconductors containing N, P, As 
and Sb in the column V of the periodic table have been extensively investigated, bismuth (Bi) containing 
alloys are less explored. Isoelectronic impurities in host materials, such as N and Bi in GaP and InP at the doping 
level, reveal rich spectroscopic features near the band gap of the host material 1 " 4 . When the impurity concen- 
tration increases from the doping level to the dilute alloy or even the concentrated alloy level, interaction between 
the impurity energy level and the band edge of the host material results in strong band gap bowing effect in dilute 
nitrides 5 ' 6 and also in dilute GaAsBi 7 . This property, which can be well explained by the band anti-crossing 
model 8 , has been employed not only to demonstrate uncooled high-speed telecom lasers on GaAs 9 and the record 
high conversion efficiency of triple junction solar cells both using dilute nitrides 10 , but also to trigger great 
intention to search for relatively unexplored dilute bismides for band engineering. Berding et al. theoretically 
analyzed InSbBi, InAsBi and InPBi as potential candidates for mid- and far-infrared (IR) optoelectronics appli- 
cations in 1988 and pointed out that InSbBi was the easiest and InPBi was the most difficult to mix, but once 
formed, InPBi was the hardest among the three materials 11 , making it an attractive narrow- gap material candid- 
ate. Since InBi crystallizes in the tetragonal PbO-type structure rather than the cubic zinc-blende structure 
common for semiconducting III-V compounds, the mixing of InBi with other materials of zinc-blende structure 
may lead to problems 12 . Early work by Oe et al. 13 and Noreika et al. 14 showed a very narrow growth window for 
InSbBi due to the weak bonding of In-Bi and formation of metallic In or Bi droplets. InAsBi samples grown at low 
temperatures show binary InBi inclusions with a tetragonal PbO phase 12 . To our knowledge, InPBi single crystal 
has never been reported since the first theoretical prediction made 25 years ago 11 . 

Here we report the first successful growth of InP x _ x Bi x thin films with Bi concentration of x = 2.4 ± 0.4% by 
gas source molecular beam epitaxy (MBE). The InPBi crystal not only reveals excellent material quality, but also 
exhibits strong and broad room temperature photoluminescence (PL) from deep levels at 1.4-2.7 um which 
provokes new theory to account for this phenomenon. 



Results 

Surface and structural properties. A set of samples with thickness of 420 nm were grown by gas source MBE and 
their growth parameters are listed in Table 1. Details of the growth procedure can be found in Methods. Sample a, 
b and c have the same flux of all three elements but different growth temperatures. Atomic force microscopy 
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Table 1 | Growth parameters and RMS roughness values mea- 
sured from AFM of the three InPBi samples grown at different 
growth temperatures 

Sample No. T g (°C) T, n (°C) T Bi (°C) PH 3 (Torr) RMS (nm) 

a 364 940 465 350 84.6 

b 324 940 465 350 0.93 

c 275 940 465 350 17.1 



(AFM) results show that droplets appear with a density in the order 
of 10 7 cm" 2 as shown in Fig. 1(a) when the growth temperature is 
364°C. The big ones are about 1200 nm in diameter and 300 nm in 
height, while the small ones are about 700 nm in diameter and 
180 nm in height. Energy dispersive X-ray (EDX) measurements 
using a scanning electron microscope (SEM) shows the relative 
percentages of In and Bi atoms to be 89.4% and 7.3%, respectively, 
in the droplets (Fig. 1(d)), while the Bi concentration is very low 
between the droplets (Fig. 1(e)). A smooth surface could be 
obtained at 324°C with a low surface roughness value of only 
0.93 nm from a 5 X 5 |im 2 scan, see Fig. 1(b). If the growth 
temperature is further lowered to 275°C, a whisker-like surface 
feature with a high roughness value of 17 nm is observed as shown 
in Fig. 1(c). The dependence of the morphology on the growth 
temperature will be discussed later. 

Figure 2(a) shows high-resolution X-ray diffraction (XRD) (004) 
00-20 rocking curves for the three samples. The right narrow peak 
with a full width at half maximum (FWHM) of around 20 arcsec 
corresponds to the InP substrate and buffer, while the left peak 
corresponds to the InPBi layer. A weak and wide peak can be seen 
in InPBi grown at 364°C. A narrow FWHM of 46 arcsec is obtained 
from the sample grown at 324°C. Interference fringes can be clearly 
seen indicating a sharp InPBi/InP interface. The InPBi thickness is 
estimated to be about 420 nm from the simulations on the fringe 
separation. The strain relaxation is determined to be less than 5.0% 
from the two symmetric (004) scans and two asymmetric (115) scans 
and the lattice mismatch with respect to the InP substrate is 0.24% 
assuming the same Poisson ratio as InP. In order to confirm the 



single crystallization of the InPBi layer, the reciprocal space maps 
(RSM) of the 430 nm thick InPBi on both (004) and (115) planes 
were performed and plotted in logarithmic scale in Fig. 2(b) and 2(c). 
Two main peaks along with several satellite peaks can be clearly seen, 
in which the one with strong intensity belongs to the substrate, and 
the one with weak intensity belongs to the InPBi layer. The FWHM of 
the InP substrate peak and the InPBi is only 12.5 and 42.0 arcsec in 
(004) map, respectively, indicating excellent single crystal quality for 
a relatively thin epi-layer. The mismatch with respect to InP is 0.24% 
with strain relaxation of 5.4% which is the same as the value obtained 
using four co-20 scans. As the growth temperature is further lowered 
to 275 °C, the InPBi peak becomes broadened. This layer has a mis- 
match of 0.25% with respect to InP and the strain relaxation increases 
to 21%. These results reveal that growth temperature is a critical 
parameter to obtain smooth InPBi single crystal. 

After the growth temperature optimized, several samples with 
similar growth parameters of sample b were grown. The Bi concen- 
tration was determined by Rutherford Backscattering Spectroscopy 
(RBS), which is a powerful method to determine the atom concen- 
tration especially for heavy atoms since the He 2+ beam will be easily 
scattered after impinging on them. A typical result is shown in Fig. 3. 
The Bi concentration is obtained by simulating the random spec- 
trum. The black dotted line, the red dashed line and the blue line are 
the randomly oriented, the simulated and the aligned spectra, 
respectively. The surface spikes of the In signal at the channel num- 
ber 455 of the [100] aligned curve is an artifact of the channeling 
analysis. Since the surface atoms are more exposed to the He 2+ beam, 
they don't channel. From the RBS results, the Bi concentration is 
determined to be 2.4 ± 0.4% with 94 ± 5% Bi atoms at substitutional 
sites. 

Figure 4(a) shows a typical cross-section Z-contrast scanning 
transmission electron microscope (STEM) image. A clear InPBi/ 
InP interface and a lightly wavy surface are observed. The InPBi layer 
thickness is measured to be 423 ± 5 nm, close to the thickness 
determined by the XRD interference fringes. From the high resolu- 
tion TEM results (not shown here), the lattice arrangement in both 
InP buffer side and InPBi side is very regular. The result indicates the 
epitaxial relation and very high interface quality. Figures 4(b) -(d) 
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Figure 1 | (a), (b) and (c) show surface morphology of 5 X 5 um 2 InPBi measured by AFM using the tapping mode when grown at 364, 324 and 275°C, 
respectively, (d) and (e) show the EDX measurement results in an SEM from the InPBi sample grown at 364°C. 
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Figure 2 | (a) High resolution XRD (004) co-29 rocking curves from the InPBi samples, (b) and (c) are the reciprocal space maps of the 430 nm thick 
InPBi on both (004) and (115) planes, respectively. 



show the EDX-mapping results in STEM for Bi, In and P, respect- 
ively. Bi incorporation in the InPBi layer is clearly revealed in quite 
uniform distribution within the probed resolution. We note, how- 
ever, information on Bi clusters or non-uniformity within the 10 nm 
regime can't be resolved from such a mapping, as the measured data 
represent the average value for each element within the e-beam probe 
which is typically of 10 nm. Nevertheless, we could not identify any 
detectable hint of Bi-related clusters from high-resolution TEM lat- 
tice imaging and high-resolution Z-contrast imaging, in contrast to 
our previous findings in GaAsBi samples 15 . For In and P mapping, 
the concentration of both elements is found to be reduced in InPBi 
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Figure 3 | Randomly oriented (black dotted line), simulated (red dashed 
line) and aligned (blue line) RBS spectra for the InPBi layer grown at 
324°C. 



compared with in InP. This is more clearly visible by the steps in 
Fig. 3(e) where depth profiles of the average Bi, In and P counts from 
EDX-mapping are shown. The detected drop of In concentration in 
InPBi is somewhat unexpected and can have important conse- 
quences for defect formations as will be discussed below. 

Optical properties. Figure 5(a) shows PL spectra of sample b at 8- 
300 K and the InP reference samples at 300 K, respectively. The 
InPBi sample at 300 K reveals strong and broad PL signals at 1.4- 
2.7 |im, centered at 1.9 (im and its FWHM is as wide as around 
660 nm (200 meV). The peak intensity is about 8 times stronger 
while the integrated PL intensity is almost 80 times higher 
compared with the InP reference sample. Careful analysis of the 
PL spectra reveals two peaks centered at 1.78 |im (0.70 eV) and 
2.04 [im (0.61 eV) at 300 K, respectively, using the Gaussian 
fitting, and three peaks at T < 200 K. This strong PL signal is not 
observed from the InP reference sample, where only an interband 
transition peak at 0.92 |im (1.34 eV) can be seen. Figure 5(b) shows 
the normalized transmission spectra of the InPBi sample and InP 
reference sample of similar thickness at 300 K. The small bump 
observed at 0.92 eV is due to water absorption. The transmittance 
drops rapidly at about 1.3 eV mainly due to the large absorption 
coefficient near the band gap of InP 16 . Compared with the InP 
reference sample, the transmittance of InPBi sample remains 
essentially the same for the photon energy below 1.0 eV, and 
decreases faster for the photon energy above 1.0 eV, indicating 
stronger optical absorption in this energy regime and a smaller 
band gap of InPBi. Using the plot of square of absorption 
coefficient vs energy (not shown here), the band gap of InPBi is 
deduced to be 1.23 eV, giving an estimate of band gap reduction of 
46 meV/%Bi. Comparing absorption with PL spectra, we notice that 
there is no energy overlaps of the interested parts of the signals. The 
absorption spectra show near band gap features at above 1.0 eV 
while the broad PL energy falls within 0.4-1.0 eV. 

Discussion 

To overcome the large miscibility gap and the weak bonding of III-Bi 
compared with those of conventional bonding of III-Vs, both the 
growth temperature and the V/Bi flux ratio should be kept low. This 
is the general strategy for almost all epitaxially grown dilute bismides. 
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Figure 4 | (a) A typical STEM image and EDX-mapping results for (b) Bi, (c) In and (d) P. Figure (e) shows the depth profiles of the average Bi, In and P 
counts from the EDX-mapping (b)-(d). 



In the case of InPBi, use of a low P flux has a risk of forming In 
droplets. For the InPBi grown at 364°C, large metallic droplets are 
observed by AFM. However, the InP thin film grown under the same 
conditions reveals a smooth surface. This indicates that the P flux is 
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Figure 5 | (a) PL spectra of the InPBi sample measured at 8-300 K and the 
InP reference sample at 300 K. The dashed blue lines are the Gaussian 
peak fitting, (b) The normalized transmission spectra of the InPBi sample 
and the InP reference sample at 300 K. 



high enough to ensure a P-rich surface reconstruction and thus a 
smooth surface. When Bi is added, large droplets are formed and the 
EDX measurement in SEM shows considerable amount of In in 
addition to Bi found in the droplets. This suggests that In is subject 
to alloying with Bi. It should be noted that both In and Bi will be in 
liquid form when the growth temperature is higher than their melt- 
ing points (156.6°C and 271°C for In and Bi, respectively) and there 
indeed exist several In-Bi phases in the phase diagram 17 . The weak 
bonding nature of Bi atoms renders themselves to accumulate on the 
InP surface. Formation of In-Bi alloys consumes both In and Bi 
atoms and hinders the InP growth. 

When the growth temperature is reduced to 324°C measured by 
the thermocouple, high quality InPBi thin film is demonstrated. 
Although the real growth temperature is unknown, we believe it will 
be quite close to the Bi melting point. The probability of Bi atoms 
sticking on surface and incorporation into lattice are enhanced when 
decreasing the growth temperature. When the impingent Bi atoms 
are balanced by their incorporation to lattice sites, single crystal 
InPBi film is achieved and no excess Bi atoms will form Bi droplets 
on surface. Bismuth is known to be as an excellent surfactant, which 
ensures good surface and interface quality. 

When the growth temperature is further lowered, excess Bi and P 
atoms will accumulate on InP. Meanwhile, PH 3 decomposition rate 
decreases with reducing growth temperature in gas source MBE and 
the resulting InPBi growth dynamics becomes rather complicated. 
The reference InP sample grown at 275°C with a slightly lower PH 3 
pressure of 320 Torr shows a similar surface feature as that of InPBi 
shown in Fig. 1(c). Further study is necessary to account for the 
observed surface morphology. 

From the (004) XRD results shown in Fig. 2(a), there is some Bi 
incorporated in the sample grown at 364°C. As the rocking-curve 
signal is broad and weak, it is difficult to estimate the Bi content. The 
sample grown at 324°C reveals nice pendellosung fringes in the (004) 
scan indicating good interface. If we assume the Vegard's law is valid, 
the extrapolated lattice constant of zinc-blende InBi will be 6.52 ± 
0.1 1 A which is very close to the previously reported value of 6.5 A 18 . 
The measured lattice mismatch is 0.24% for InPBi grown at 324°C. 

The PL measurement reveals unexpectedly strong and broad sig- 
nals at 0.4-1.0 eV at 300 K which are absent in the reference InP thin 
film grown at the same growth conditions. The absorption spectra 
show no difference between InPBi and InP for energy below 1.0 eV 
but enhanced absorption in InPBi for energy above 1.0 eV. This 



SCIENTIFIC REPORTS | 4:5449 | DOI: 1 0.1 038/srep05449 



4 



implies that physical origins can be different in the two energy 
regimes. Previous spectroscopic studies of InP:Bi grown by liquid 
phase epitaxy with the Bi concentration in the doping level of 10- 
100 ppm show many spectral features within 100 meV below the 
InP band gap at low temperatures, but no signals in the mid-gap at 
room temperature 3,4 . In these materials, Bi can be regarded as 
dopant and modelled by a resonant defect level within the energy 
band structure of the host matrix. In our case, the Bi concentration 
is much higher, and Bi interaction as well as cluster formation 
must be considered. There have been several theoretical investi- 
gations on the band structures of GaAsBi 19-23 and GaPBi 23 to account 
for the band gap reduction of dilute bismides, but the physical 
explanation of whether it is due to the band anti- crossing or just 
the alloying and disorder is still under debate 23 . It is theoretically 
predicted that the Bi-related defect resonant level lies 80 meV 
below the valence band of GaAs 22 and 122 meV above the valence 
band of GaP 23 . For disordered GaP x _ x Bi x including clustering 
effect with x up to 4.88%, broad band tail states exist at about 
0.15-0.5 eV above the valence band 23 . If we assume similar band tail 
states also exist in InPBi, this may explain the near band gap features 
observed in the absorption measurement. However, it gives low 
energy bound of 0.84 eV (the band gap of InP is 1.34 eV at 
300 K), still larger than the observed 0.70 eV and 0.61 eV in PL at 
300 K. On the other hand, if the PL signals were due to the band anti- 
crossing, this would give exceptionally large band gap bowing of 
267 meV/%Bi. Therefore, current existing theory can't explain the 
PL signals at the mid band gap range. Reihlen et al. reported large 
energy differences of 70-200 meV between PL and absorption at 
300 K in InP x Sb! _ x for O.K x < 0.4 24 . They attributed this differ- 
ence to the combination of deep centers and intrinsic band tail 
recombination which were indistinguishable from the spectroscopy. 
The large band tail was believed to result from the compositional 
clustering. The compositional inhomogeneity usually induces loca- 
lized states and causes S-shape in temperature dependent PL peak 
energy, often found in dilute nitrides 25 . However, in the temperature 
dependent PL spectra shown in Fig. 4(a), we don't observe any 
noticeable S-shape curves. 

We tentatively attribute the strong and broad PL signals at 300 K 
to deep level luminescence. The deep levels can be Bi In anti-sites, P 
and/or In vacancies or the above combinations. The EDX- mapping 
results shown in Fig. 3 unambiguously confirm the decreases of both 
In and P atoms in InPBi compared with in InP. This gives some 
evidence that In and P vacancies may exist. Indeed, the scanning 
tunneling microscope (STM) on cross -sectional InPBi confirms the 
existence of P vacancies (not shown here). Bi atoms are generally 
believed to occupy P sites in dilute bismide. However, this will intro- 
duce very high local strain in InPBi. Our first principle calculations 
show that the overall formation energy is smaller for Bi taking In 
lattice positions instead of P positions. Also Bi can be of 3 + valence 
due to the relativistic effect of the two 6 s electrons. These suggest 
that Bi In anti- site is physically favorable and its concentration can be 
quite high. More experimental works, such as deep level transient 
spectroscopy (DLTS) and time resolved PL, are necessary to clarify 
the nature of those deep levels and carrier dynamics. Likewise, new 
theoretical model or simulations are stimulated to give physical 
explanations of the observed PL spectra. 

In summary, new InPBi single crystal has been successfully grown 
by gas source MBE for the first time. High surface and structural 
quality has been demonstrated and confirmed by AFM, XRD and 
TEM measurement. RBS measurement shows that the Bi concentra- 
tion is 2.4 ± 0.4%, with 94 ± 5% Bi atoms at substitutional sites. 
Room temperature PL measurement reveals a very broad and strong 
spectrum from 1.4 to 2.7 um while the optical absorption shows a 
band gap reduction of 1 10 meV with respect to InP. The observed PL 
spectrum is likely from deep levels and can't be explained by the 
existing theory. 



Methods 

Molecular Beam Epitaxy. Samples were grown on (100) semi- insulating InP 
substrates by V90 gas source MBE. P 2 was cracked from phosphine at 1000°C and its 
flux was controlled by regulating the gas pressure in the gas line. Elementary In and Bi 
sources were used and their fluxes were controlled by adjusting the effusion cell 
temperatures. Both substrate and cell temperature were measured by thermocouples. 

Prior to growth, oxide desorption from the substrate surface was carried out by 
heating the substrate with a P 2 flux impingement on the surface. The complete 
desorption of the oxide layer was determined using the reflection high energy electron 
diffraction (RHEED), by a sharp transition to (2 X 4) pattern at about 524°C. A 
75 nm un-doped InP buffer layer was grown at 474°C with the PH 3 pressure of 630 
Torr to improve the surface quality. 420 nm thick InPx _ x Bi x thin films were then 
grown at 275-364°C at 0.86 um/h with a low P flux to facilitate the Bi incorporation 
into the lattice. InP reference samples were grown with otherwise the same growth 
conditions as those for InP : _ x Bi x for comparison. 

Measurements. The surface morphology was measured by Bruker Icon atomic force 
microscopy (AFM) using the tapping mode. The scan area was 5X5 um 2 . The 
electron acceleration voltage for the EDX measurement in SEM was set at 10 kV to 
obtain good spatial resolution. Conventional X-ray diffraction co-20 scan was 
performed using a Philips X-pert high -resolution diffractometer to measure 
structural properties. The Cu KM wavelength was selected by a Ge (220) fourfold 
monochrometer. The STEM and the EDX- mapping measurements were performed 
on a FEI Tecnai G2 F20 microscope. Cross- section STEM specimens were prepared 
by focus ions beam. In the RBS measurement, a 2.275 MeV 4 He 2+ source was used to 
impinge on the sample surface. The PL and transmission spectra were measured using 
a Nicolet Megna 860 Fourier transform infrared (FTIR) spectrometer, in which a 
liquid- nitrogen cooled InSb detector and a CaF 2 beam splitter were used. A diode- 
pumped solid state laser (X = 514.5 nm) was used as an excitation source for PL 
measurement. 
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